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THE LIFETIMES OF MOLECULAR CLOUD CORES: WHAT IS THE ROLE 

OF THE MAGNETIC FIELD? 
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RESUMEN 

Discutimos los tiempos de vida y la evolution de los nucleos densos formados como fluctuaciones turbulentas 
de densidad en nubes moleculares isotermicas magnetizadas. Consideramos simulaciones numericas en las que 
mcdimos la criticalidad magnetica y la estabilidad de Jeans de los nucleos, en relation a la criticalidad magnetica 
de sus nubes "madre" (los dominios de integration numcrica). En dominios subcriticos, no sc forman nucleos 
densos, y no ocurre colapso gravitacional. En dominios supercriticos, algunos nucleos se colapsan, formando 
parte de grumos mas grandes que son supercriticos desde el initio, y cuyas regiones mas densas (los nucleos) 
son inicialmcntc subcriticas, pero rapidamcnte se tornan supercriticos, presumiblcmente por acrecion a lo largo 
de las lineas de campo magnetico. Dcscartamos la posibilidad de artefactos numcricos. Las cscalas de tiempo 
en las que los nucleos transitan desde el cstado subcritico, pasando por el supcrcritico, hasta finalmente el 
colapso, son de unas cuantas veces su tiempo de caida librc, rf c . Nucstros rcsultados sugieren que los nucleos 
son estructuras transientes, fuera de equilibrio, y no configuraciones cuasi-magneto-hidrostaticas. 

ABSTRACT 

We discuss the lifetimes and evolution of dense cores formed as turbulent density fluctuations in magnetized, 
isothermal molecular clouds. We consider numerical simulations in which we measure the cores' magnetic 
criticality and Jeans stability in relation to the magnetic criticality of their "parent clouds" (the numerical 
boxes). In subcritical boxes, dense cores do not form, and collapse does not occur. In supercritical boxes, some 
cores collapse, being part of larger clumps that are supercritical from the start, and whose central, densest 
regions (the cores) are initially subcritical, but rapidly become supercritical, presumably by accretion along 
field lines. Numerical artifacts are ruled out. The time scales for cores to go from subcritical to supercritical and 
then collapse are a few times the core free-fall time, Tf c . Our results suggest that cores are out-of-equilibrium, 
transient structures, rather than quasi-magnetostatic configurations. 
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1. INTRODUCTION 

The prevailing view (which we hereafter refer to 
as the "standard (magnetic support) model" of star 
formation; see, e.g., the reviews by Shu, Adams & 
Lizano 1987; McKee ct al. 1993) concerning low- 
mass-star-forming clumps is that they are quasi- 
static equilibrium configurations with so-called "sub- 
critical" mass-to-magnctic-flux ratios, so that the 
clumps are supported against their self-gravity by 
the magnetic field in the direction perpendicular to 
it, and by a combination of thermal and turbulent 
pressures along the field. Under ideal MHD condi- 
tions, the magnetic field is "frozen" into the plasma, 
and the magnetic flux is conserved. Under the addi- 
tional assumption that the clump's mass is also con- 
stant, then the mass-to-flux ratio is a fixed paramc- 
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ter of the clump, which therefore cannot collapse if 
this ratio is subcritical (meaning that the core's self- 
gravity is never enough to overwhelm the magnetic 
support). However, because the cold molecular gas 
is only partially ionized, the process known as am- 
bipolar diffusion causes a loss of magnetic flux from 
the clumps on time scales long compared to their 
free-fall time, allowing them to contract and form 
denser cores that will ultimately collapse. 

However, it well known that molecular clouds are 
supcrsonically turbulent (e.g., Larson 1981; Blitz & 
Williams 1999), and it is becoming increasingly ac- 
cepted that the cores within them are the density 
fluctuations induced by the turbulence (Ballesteros- 
Paredes, Vazquez-Semadeni & Scalo 1999; see also 
the reviews by Vazquez-Semadeni et al. 2000; Mac 
Low & Klessen 2004). In this context, the cores have 
a highly dynamical origin (supersonic compressions), 
and their masses are hardly fixed. Moreover, it is 
natural to ask whether they can settle into hydro- 
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static equilibria, an event which requires the equi- 
libria to be stable (or "attracting" , in the language 
of nonlinear phenomena). Otherwise, the dynamic 
density fluctuations will just "fly past" the equilib- 
rium state on their way to collapse, or else "rebound" 
and merge back with their environment, if they do 
not quite reach the equilibrium point. In this case, 
the cores' lifetimes should be much shorter than in 
the standard model, probably comparable to their 
free-fall times. In the present contribution, we ar- 
gue in favor of this scenario. To this end, we discuss 
the formation and evolutionary time scales of cores 
that form in numerical simulations of isothermal, 
compressible MHD turbulence, in relation to the 
magnetic criticality of the whole computational box. 
Here we present a brief overview. For a full, detailed 
discussion, see Vazquez-Semadcni ct al. (2004). 

2. QUALITATIVE CONSIDERATIONS 
For the discussion below, it is convenient to char- 
acterize a core by two non-dimensional parameters: 
its Jeans number J c = R c /Lj, giving the ratio of the 
core's radius R c to its Jeans length ij, so that a 
gravitationally unstable core (with respect to ther- 
mal support) has J c > 1, and its mass-to- flux ratio 
/i c , normalized to the critical value for magnetic sup- 
port, so that a supercritical structure has fi c > 1. 
The corresponding values for the entire computa- 
tional box are denoted J and /i. 

To assess the necessary conditions for the for- 
mation of either sub- or supercritical cores, we note 
that the mass-to-flux ratio of an isolated cloud con- 
stitutes an upper bound for that of any subregion 
within it, since, as long as the flux-freezing condi- 
tion holds, both the cloud's mass and magnetic flux 
are fixed, regardless of the cloud's volume. So, even 
if the entire mass of the cloud were compresed into 
the core's volume, its mass-to-flux ratio would re- 
main the same. Thus, in the absence of ambipolar 
diffusion, supercritical cores can only form within su- 
percritical clouds or clumps. On the other hand, sub- 
critical cores can arise in either sub- or supercritical 
clouds. Now, a gravitationally bound, yet subcritical 
clump must have J c > 1 in addition to fx c < 1, so 
that it does not re-expand after the turbulent com- 
pression that formed it subsides. However, recent 
observational studies of the magnetic field strengths 
in molecular clouds suggest that these are the ob- 
jects with highest mass-to-flux ratio in the hierar- 
chy going from diffuse clouds to dense cores (e.g., 
Crutchcr 2004), and are likely to be supercritical in 
general (e.g., Bourke et al. 2001), so the formation 
of supercritical cores without the need for ambipolar 
diffusion appears feasible in real molecular clouds. 




Fig. 1. Evolution of the global maximum of the density 
field for all runs considered here. The numbers indicate 
the corresponding values of (3. 

3. NUMERICAL SIMULATIONS AND RESULTS 

In order to investigate the magnetic criticality of 
the cores that can form in turbulent environments in 
relation to that of their parent clouds, we have per- 
formed a suite of numerical simulations of turbulent 
flows at fixed rms sonic Mach number (M s ~ 10) and 
global Jeans number J = 4, and varying the plasma 
(3, defined as the ratio of thermal to magnetic pres- 
sure, in order to consider cases that are subcritical 
(denoted /3.01, with fx = 0.9), mildly supercritical 
(denoted (3.1, with fi — 2.8), strongly supercritical 
(denoted (31, with /i = 8.8), and non-magnetic (de- 
noted Poo, with fi — oo). These simulations can be 
thought of as representing regions of size L = 4 pc, 
with a mean number density no = 500 cm~ 3 , a tur- 
bulent velocity dispersion of 2 km s _1 , and a sound 
speed c s = 0.2 km s _1 , with mean field strengths of 
46, 14.5, 4.5 and uG, respectively. 

Figurcnshows the evolution of the global density 
maximum for all four runs, with the time axis shown 
in units of the sound crossing time r s = L/c s = 20 
Myr (lower axis) and of the global free-fall time 
Tf g = Lj/cg = 5 Myr (upper axis). It is readily seen 
that the subcritical run (3.01 does not produce very 
large density enhancements, with the global density 
maximum at any given time being ~ 30no, and never 
exceeding 100?io- These values are too low, and the 
transients are too short (< 1 Myr), for ambipolar 
diffusion to operate and reduce the magnetic flux 
under canonical estimates of the ambipolar diffu- 
sion time scale (sec, e.g., McKcc ct al. 1993). Note 
also that this occurs even though this run is very 
close to global criticality, and thus the production 
of gravitationally-bound cores in subcritical environ- 
ments appears to be highly unlikely even for nearly 
critical boxes. Finally, note that the densities en- 
countered in this run are also safely below the "Jeans 
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criterion" of True-love et al. (1997) (n < 256no for our 
resolution of 256 3 grid zones) and the extension of 
it proposed by Heitsch, Mac Low & Klessen (2001) 
(n < 115no) to ensure that numerical diffusion does 
not significantly damp MHD waves within the cores, 
and so our result is robust. Note, moreover, that 
these criteria may be excessive for the problem at 
hand, as it was shown by Heitsch et al. (2001) that 
the magnetic field can only prevent collapse if it 
provides magneto static support, rather than wave- 
pressure support. So our main resolution concern is 
to avoid magnetic flux loss, not damping of MHD 
waves. In turn, the Truelove et al. (1997) criterion 
is only defined in order to prevent artificial fragmen- 
tation of collapsing structures, but here we are not 
concerned with the internal structure of the cores 
as they collapse; rather, we are just interested in 
whether collapse occurs or not. Thus, a more appro- 
priate criterion may be that the internal mass-to-flux 
ratio of the cores does not exceed that of the numeri- 
cal box, or of its parent structure, as discussed above 
(i.e., that the mass-to- flux ratio does not increase in- 
wards of a density structure). 

In contrast, the supercritical and nonmagnet- 
ic runs rapidly develop densities approaching 
10 no, which correspond to collapsed, unresolved 
objects. These densities occur at ~ 1/2t{ & in 
the magnetic runs, and at ~ l/5rf g in the non- 
magnetic one. However, inspection of animations 
of the simulations (available at http://www. 
astrosmo . unam . mx/ ~e . vazquez/turbulence_HP/ 
movies/VKSB04 . html, or in Vazquez-Semadcni et 
al. 2004) shows that the typical time spans from 
the beginning of the turbulent compression to the 
completion of collapse are ~ 1-1.5 local (i.e., at the 
mean density of the core) free-fall times (Tf c ), or, 
equivalently, ~ 0.1-0.15 Tf g . 

Measurement of the Jeans number and mass-to- 
flux ratio of the cores as they are formed and evolve 
towards collapse sheds light on how the process oc- 
curs. In Table ^ we present these parameters, to- 
gether with other relevant data, for the first collapsed 
object that forms in run /3.1. Other cores have sim- 
ilar histories. The animation shows that this object 
forms out of a larger clump initially containing two 
cores, which ultimately merge to form the collapsed 
object (see also Figure |5J). Table Q] thus gives the 
data for the structures defined out to a threshold 
density level n t from the local maxima at three dif- 
ferent times during the evolution of the system. The 
times are in units of the box sound-crossing time. At 
the time t = 0.06t s , the threshold n t = 40no resolves 
the two cores, but at the two later ones only one core 




Fig. 2. Iso-density surface 3D map of the clump-core 
(cores 1 and 2) system appearing at early times in run 
/3.1, detailed in Tabled The times shown are t — 0.06r s 
(left) and t = 0.08r s (right). The iso-density surfaces are 
at n = 10no (blue), n = 40no (yellow) and n = 100no 
(red). The latter does not exist at t = 0.06r s . 

is seen above this threshold (see Figure 

From Table Q it is seen that the parent clump 
(defined by setting n t = 10uq) is already super- Jeans 
and supercritical at t = 0.06r s , but the daughter 
cores (defined by n t = 40no) are still sub-Jeans and 
subcritical. Both the mean and peak densities (re- 
spectively n and n p ) within the parent and daughter 
structures are safely below even the most stringent 
condition of n < 115no for avoidance of MHD wave 
damping by numerical diffusion, and so the super- 
critical nature of the parent clump is a robust physi- 
cal result, and the clump cannot be supported by the 
magnetic field. Of course, there remains the possibil- 
ity of turbulent support, but since the clump itself 
was formed by a turbulent compression, the com- 
pressive energy clearly overwhelms the internal, ran- 
dom, supporting one, and in fact is the mechanism 
that drives the clump into becoming super- Jeans and 
supercritical. In any case, a detailed analysis of the 
virial balance of the clump/core system and the role 
of the velocity field, to be presented elsewhere, is 
necessary, but here we have shown that the magnetic 
field is insufficient for supporting the parent clump, 
even though its densest regions (the daughter cores) 
are initially subcritical and sub-Jeans. At the later 
times t = 0.08t s and t = 0.1t s , the single core defined 
by n t = 40 also has become super- Jeans and super- 
critical, as expected for the generalized collapse of 
the parent clump. Although finally at time t = 0.1t s 
numerical diffusion is clearly important (fi c > u), 
this is seen to occur after the onset of the clump's 
collapse, which began while the core/clump system 
was still well resolved. 

A final observation that here we only mention 
briefly (see Vazquez-Semadeni et al. 2004 for a full 
discussion) is that in the strongly supercritical run 
/31, two cores with slightly longer durations (~ 5Tf c ) 
are formed, but in both cases they end up re- 
dispersing, rather than collapsing. This is in agree- 
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TABLE 1 

PARAMETERS OF CLUMPS AND CORES IN RUN [3.1. 



Core name 


Time a 


n t b 


n c 


n p d 


R c /L c 


/ f 


A* c s 


Parent clump 


0.06 


10 


19.0 


83.0 


0.071 


1.27 


1.73 


Core 1 


0.06 


40 


51.1 


83.0 


0.018 


0.53 


0.90 


Core 2 


0.06 


40 


44.4 


56.0 


0.020 


0.54 


0.86 


Single core 


0.08 


40 


78.5 


337. 


0.041 


1.45 


2.09 


Single core 


0.10 


40 


329. 


5.07 x 10 3 


0.031 


2.28 


4.66 



a Timc of observation of structure, in units of t s . 
b Threshold density for clump/core identification. 
c Mean density of clump/core. 
d Maximum density within clump/core. 

c Clump/core radius in units of the computational box size L. 
f Ratio of core's radius to its Jeans length. 

s Mean mass-to- flux ratio within clump/core in units of the critical value (compare to simulation value of 



ment with the simple argument that re-expansion 
should take longer times than compression or col- 
lapse, because in that case self-gravity acts as re- 
tarding agent against the re-expansion. 

4. CONCLUSIONS 

In this contribution, we have presented a dis- 
cussion and numerical simulations of the formation, 
nature, and lifetimes of dense cores in magnetized 
clouds. We argued that the mass-to-flux ratio of a 
cloud puts an upper bound to that of any clump 
or core within it, and so the only way to form su- 
percritical cores is within supercritical clouds (ne- 
glecting ambipolar diffusion). Moreover, numeri- 
cal simulations of marginally subcritical clouds show 
that no gravitationally bound cores form in this 
case, while in supercritical clouds the gravitation- 
ally bound cores that form occur inside clumps that 
are supercritical, and rapidly become supercritical 
themselves and collapse, in timescales that do not 
exceed twice their local free-fall time Tf c . We also 
showed that the mass-to-flux ratio of structures de- 
fined through a density threshold (as would be the 
case of cores observed in a single molecular line, 
which requires the density to be larger than a certain 
value in order to be excited) does not remain con- 
stant, because the core is continuously connected to 
its parent structure, from which it can accrete mass. 
Our results support the notion that clumps and cores 
are out-of-equilibrium, transient structures, and that 
a class of "failed" cores should exist, that will not 
form stars. 
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